Dephasing due to electron-electron interaction in a diffusive ring 
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We study the effect of the electron-electron interaction on the weak localization correction of a ring 
pierced by a magnetic flux. We compute exactly the path integral giving the magnetoconductivity 
for an isolated ring. The results are interpreted in a time representation. This allows to characterize 
the nature of the phase coherence relaxation in the ring. The nature of the relaxation depends on the 
time regime (diffusive or ergodic) but also on the harmonics n of the magnetoconductivity. Whereas 
phase coherence relaxation is non exponential for the harmonic n = 0, it is always exponential for 
harmonics n ^ 0. Then we consider the case of a ring connected to reservoirs and discuss the effect 
of connecting wires. We recover the behaviour of the harmonics predicted recently by Ludwig & 
Mirlin [Phys. Rev. B69 (2004) f 93306] for a large perimeter (compared to the Nyquist length). We 
also predict a new behaviour when the Nyquist length exceeds the perimeter. 

PACS numbers: 73.23.-b ; 73.20.Fz ; 72.15.Rn 



I. INTRODUCTION 

In the classical description of transport in weakly disor- 
dered metals, elastic scattering by impurities leads to the 
finite Drude conductivity at low temperature. It is well- 
known that quantum interferences manifest themselves 
through a small sample dependent contribution, whose 
average, denoted by (Act), is called the weak localization 
correction. Dephasing strongly affects weak localization, 
which provides a powerful tool to probe phase coherence 
in disordered metals. The simplest approach to describe 
dephasing is to assume that the time dependence of the 
phase coherence relaxation is exponential. Such relax- 
ation can be due for example to magnetic impuritiesi. It 
is characterized by a time scale called the phase coher- 
ence time T v and the weak localization correction to the 
conductivity takes the form 
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where the so-called cooperon V(t) is the contribution to 
the return probability originating from quantum inter- 
ferences between time reversed trajectories. It is solu- 
tion of a diffusion equation. The factor 2 stands for spin 
degeneracy and D is the diffusion coefficient. We have 
set h = 1. In a quasi-one-dimensional infinite wire, the 
probability is well known to vary as V(t) — 1/ S\J AnDt, 
where S is the cross section of the wire, so that the weak 
localization correction has the familiar form 
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where L v = \jDi\p is called the phase coherence length. 

The measurement of the weak localization correction 
is possible thanks to its sensitivity to an external mag- 
netic field. For a wire^, the effect of a weak perpendic- 
ular magnetic field can be described by introducing an 



exponential reduction factor e 



-t/T B 



JU , where the 



m eq. 

characteristic time is tb = 3/(De 2 B 2 S) (for a wire of 
square cross section). Consequently the weak localiza- 
tion is given by eq. J2J with the addition of the inverse 
times, "a la Matthiesen" : 
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From the experimental point of view, this effect is of 
primary importance, since the magnetic field acts as a 
probe in order to study phase coherence and to extract 
t v and its temperature dependence. 

In the more complicated geometry of a ring, the field is 
also responsible for magnetoconductivity oscillations as 
predicted by Altshuler, Aronov & Spivak (AAS) 3 . The 
phase coherent return probability is sensitive to the flux (f> 
through the ring. It has the simple harmonics expansion : 
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where 9 = 47r</>/(/>o is the reduced flux (</>q = h/e is the 
flux quantum). Each harmonic corresponds to a number 
of windings of the diffusive trajectories around the ring. 
The relation (JTJ, with 10}, immediately leads to the fa- 
miliar result of AAS for the weak localization correction 
to the average conductivity in a ring : 



(Aa(9)) = 



e 2 sinh(£/L y ) 
irS v cosh(L/L v ) — cosO 



(5) 



The harmonics of this ^o/2-periodic correction decay ex- 
ponentially with the perimeter L of the ring : 

(Act„) = / — (Aa(9)) e~ ine = L v e-l"IVV(6) 

The combination of the two effects of the magnetic field, 
AAS oscillations and penetration in the wires, is obtained 
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by perfoming the substitution : 1/t v — > 1/t v + 1/tb in 
eq. ©. 

Despite the exponential damping in eq. (JIJ describes 
correctly several dephasing mechanisms like spin-orbit 
scattering and spin-flipi, or the effect of an external 
magnetic field^, a precise description of the electron- 
electron interaction requires a more elaborate treatment. 
In a pioneering paper, Altshuler, Aronov & Khmelnit- 
skii (AAK) 4 *^ have shown that the dephasing due to the 
electron-electron interaction can be described in a one- 
particle picture by coupling the electron to a fluctuating 
classical electromagnetic field. They obtained a result 
which can be cast in the form 

(Act) = -2 / dtV(t) f(t/r N ) e" 4 ^ , (7) 
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where f(x) is a decreasing dimensionless function. We 
have also incorporated an exponential relaxation (if it 
is due to the external magnetic field B we have simply 
Tip — > tb). The characteristic time, called the Nyquist 
time, is given by& 
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where T is the temperature and ks the Boltzmann's 
constant (in the following we will set H — ks = 1). 
(To = 2e 2 pgD is the Drude conductivity and po the den- 
sity of states at Fermi energy, per spin channel. The 
y-2/ 3 power law has been observed in a variety of ex- 
periments (see for example refsAS) and is the signature 
of the electron-electron interaction in quasi-ld wires. In 
the case of an infinite wire, AAK found that the weak 
localization correction is given by^ 



(Act) = —L 



Ai(t n /t v ) 

nS Ai' (t n /t v ) 
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where Ai(z) is the Airy function and Ai (z) its derivative. 
We have introduced the Nyquist length Ljy = VDt~j 
which characterizes the scale over which the electron- 
electron interaction is effective. It can be conveniently 
rewritten as : 
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where we have introduced the thermal length Lt — 
y/D/T, the number of channels N c , the elastic mean 
free path £ e and a numerical factor aj, that depends 
on the dimension d (a% — 2, ai = n/2 and — 
4/3) 10 i 11 ' 12 . We have expressed the Drude conductiv- 
ity as cjq = 2^-adN c £ e / S (the factor 2 stands for spin 
degeneracy) . 

In addition to the prediction of the power law tn cx 
T~ 2 / 3 for the coherence time, an important outcome of 
the AAK theory is that the result JjJJ obviously breaks 



the addition rule of inverse times. This indicates that the 
phase relaxation characterized by the function /(i/rjv) 
is indeed non exponential. This function has been calcu- 
lated recently in refi^ where it was found that it varies 



14 p-( % A?/4)( t /r JV ) 3/2 



for t < Tjv- However the study 



as 

of the function /(£/ T iv) shows that it is very close to 
an exponential e~ t / 2TN and eq. © only deviates from 

(Act) = -^#(5^7 + ^)- 1/2 , that is eq. ©, by no 
more than 4%i2 i i&± 6 .. This explains why it is very dif- 
ficult to observe experimentally the functional form 
and most of the magnetoconductance measurements in 
wires have been analyzed assuming the form ©. 

In the paper of AAK, only simple geometries (like wire 
and plane) were considered and it is not clear how the 
non exponential relaxation of phase coherence affects the 
weak localization for a nontrivial geometry. In a recent 
paper, Ludwig & Mirlin (LM)i£ have addressed the ques- 
tion of dephasing due to the electron-electron interaction 
in a ring. The dephasing is then probed by the harmonics 
of the magnetoconductance oscillations. These authors 
found that these harmonics decay with the perimeter L 
of the ring in an unexpected way. LM's result can be cast 
in the formic 
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This result is quite interesting, because this nontrivial 
non exponential decay of the harmonics leads to an un- 
expected e _ni3/2Tl/2 temperature behaviour, instead of 
the incorrect behaviour e _ni / ijv = e~ nLT /3 naively ex- 
pected from a simple substitution t v — > r^r in the AAS 
harmonics iJBJ. It also shows that the geometry of the 
system may play an important role in the nature of the 
dephasing mechanism. 

The work of LM was mainly devoted to the study of 
Aharonov-Bohm (AB) oscillations in a single ring. The 
study of AB amplitude, rather than AAS, is motivated 
by the lack of disorder averaging in this case. The ampli- 
tude of AB oscillations is given by the harmonics (Sa^) of 
the conductivity correlation function (8a(B) 6<r(B')). As 
pointed out by LM, these harmonics are expected to be 
directly related to the AAS harmonics by the following 
relation 



(So* 



(Act* 
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where Lt = \J D/T is the thermal length. This expres- 
sion extends the result of Aleiner & Blanter— who stud- 
ied the relation between conductance fluctuations and 
weak localization in a wire and a plane when dephasing 
is due to the electron-electron interaction. An important 
consequence of this relation is that the effect of dephas- 
ing on weak localization and conductance fluctuations is 
governed by the same length scale L^- We re-examine 
this relation and give a more general proof in appendix lEl 

In our paper, we reconsider the question of weak local- 
ization in a ring in the presence of electron-electron inter- 



3 



action. Our main goal is to provide a physical picture as 
well as a detailed understanding of the results obtained 
by LMii. The physical reason for the geometry depen- 
dence of the dephasing can be understood in the following 
heuristic way. For a pair of time reversed trajectories, we 
denote by <& the random phase brought by the fluctuating 
electromagnetic field. Average over the Gaussian fluctu- 
ations of the field is denoted by (• • -)v- Averaging the 
phase (e 1 *)^ produces an exponential damping e - '^* ^ v 
responsible for phase coherence relaxation (this exponen- 
tial is related to the function /(t/rjv) m eq. Q0- F° r a 
quasi-ld system, the typical damping rate associated to 
a diffusive trajectory can be written in the form 



d($ 2 



dt 
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where r(t) designates the typical distance explored by 
the diffusive trajectory over a time scale t. As pointed 
out in refi^S, eq. lfl3"j) can be understood as a local form 
of the Johnson-Nyquist theorem : since the phase and 
the potential are related by $ = V, eq. (|T3|) measures 
the potential fluctuations ^($ 2 )v = f dt(V(t)V{0)) v = 
2e 2 TR t , where Rt ~ r(t)/(S<Jo) is the resistance of a wire 
of length r(t). It is clear from eq. I|13(l that the dephasing 
depends on the nature of the diffusive trajectories. Two 
regimes can be distinguished : 

• The diffusive regime.- In this regime, the bound- 
aries of the system play no role and the diffusion 
follows the behaviour obtained in an infinite wire, 
therefore r(t) ~ VTJt, so that 



3/2 



(14) 



Therefore we expect the phase relaxation to be non 
exponential /(t/Tjv) ~ exp — (t/"nv) 3 / 2 . 

• The ergodic regime. - In this case, the diffusive tra- 
jectories explore the whole system. This corre- 
sponds to time scales larger than the Thouless time 
to = L 2 /D. The characteristic length is given by 
the size of the system (the ring) r(f) ~ L and 



/^2\ V^D^ * 
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We expect the phase relaxation to be exponential, 
exp— t/r c , where the time scale 
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is size dependent and is a nontrivial combination of 
the Thouless time and the Nyquist time. 

It is clear that the function /(t/rjv) in eq. Q must be re- 
placed by a more complicated function f (t / tn , t / tjj) to 



account for finite size effect and describe both regimes. 
We expect that the winding around the ring plays an 
important role and a dephasing of different nature for 
trajectories which enclose the ring (corresponding to the 
n 7^ harmonics of the flux dependence) and for trajecto- 
ries which do not encircle the ring (corresponding to the 
n = harmonic). For harmonics n ^ 0, the trajectories 
necessarily explore the whole ring and the phase coher- 
ence relaxation is always exponential. This is the origin 
for the new behaviour found by LM. For the harmonic 
n = we can distinguish two different regimes corre- 
sponding to non exponential and exponential relaxation. 
In other terms, the function /(t/r^r,t/rc) depends also 
on the winding number n. 

In order to give a firm basis to these arguments, we 
have reconsidered the calculation of LM. These authors 
have studied a symmetric ring connected by two arms to 
reservoirs and have derived the weak localization correc- 
tion within an instanton approximation for the functional 
integral describing the effect of the fluctuating field. The 
instanton approach is only valid in the limit of large 
perimeter of the ring (compared to the length scale char- 
acterizing the electron-electron interaction), that is when 
L > L N {i.e. t d > t n ). 

In our work we have followed a different strategy. We 
first consider the case of an isolated ring for which it 
is possible to compute exactly the path integral for any 
value of L/Ln. This path integral formulation is intro- 
duced in section [HI and the exact solution is given in 
section ITTTl where a closed expression for the harmonics 
of the magnetoconductivity is provided and analyzed in 
several regimes. 

Our result agrees with the exponential behaviour i|ll|) 
found by LM. However we will point out that LM es- 
timated an incorrect prefactor (leading to an incorrect 
temperature dependence) . LM's result can be written as 
(Aer n ) LM ~ L^ 4 e" |n|(i/ - Ljv:)3/2 whereas we will show that 

the correct result is (Aa n ) ~ J L A re _|n|(L/Ljv)3/2 . Note 
that it could seem at first sight that the difference comes 
from the fact that our exact solution stands for an iso- 
lated ring whereas LM considered a connected ring, how- 
ever we will see that the two situations are closely related. 
Moreover the isolated ring leads to consider a path in- 
tegral exactly similar to the one studied by LM which 
allows us to trace back to the origin of LM's incorrect 
prefactor (we have also estimated in appendix[C]the path 
integral within the instanton approach used by LM. We 
treat carefully the prefactor within this approach). 

Our exact solution takes also into account the effect of 
an exponential relaxation of phase coherence (t v ). Wc 
will show that both kinds of dephasing mechanisms (i.e. 
tjv and t v ) combine in a nontrivial way. 

In section Hvl we analyze the exact result in a time rep- 
resentation, that is we study the function /(t/rjvj £/td) 
that generalizes the /(i/ 7 "^) 01 ec l- 0- A special em- 
phasis is put on the difference between harmonics n = 
and n^O : although the phase coherence relaxation for 
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the harmonic n = is either non exponential (at short 
times) or exponential (at large times), it is always ex- 
ponential for harmonics n =/= 0. This analysis in time 
representation is used in the section Ivl where we discuss 
the effect of connecting wires : in a transport experiment, 
the ring is necessary connected to wires that can strongly 
affect the magnetoconductance. This has been discussed 
recently for the case of exponential relaxation of phase 
coherence^. Although we do not expect a strong effect 
of the connecting wires on the harmonics in the regime 
L Ln, we will show by some simple arguments that the 
behaviour is strongly modified in the other regime 
L < L N . 



II. PATH INTEGRAL FORMULATION 

We recall the basic ideas of AAK's approach^. In di- 
mension d ^ 2 the dephasing is dominated by small en- 
ergy transfers. The dephasing for one electron can be 
modeled through the coupling of the electron with a clas- 
sical fluctuating electric potential V(r,t), whose fluctu- 
ations are given by the fluctuation-dissipation theorem. 
In a Fourier representation, the correlations are given 
by (VV)(q^) — l^pr- In a time-space representation it 
reads 



(V(f,t)V(f\l/)) v = 2 -^T8{t-t')P d {r,r') , (17) 



where 



Co 



designates averaging over the Gaussian fluc- 



tuations of the potential. Pd{r,r') is the diffuson, solu- 
tion of the equation —APd(r, r') = S(r — r') (it is under- 
stood that the zero mode diverging contribution is not 
taken into account in the diffuson of eq. \VT\- This point 
is discussed in appendix^). Eq. (fP7|) clearly exhibits the 
long range nature of the spatial correlations. 

The description of AAK allows to perform a mean- 
field-like but non-perturbative treatment of the electron- 
electron interaction. We denote by V the cooperon that 
includes the effect of the electron-electron interaction. It 
can be conveniently written with a path integration over 
the diffusive trajectories of the electron 



^r{i)—r 



T>t(t) e 
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+ i*[r(r)] 
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where the phase 

$[f%r)] = ( dr \V{?{t),t) - V(r(r),f)] , (19) 
Jo 

with t — t — t, is the phase difference between the two 
time reversed trajectories. In the presence of a static 
magnetic field, a coupling 2ie L drf^r) • A(r(r)) to the 
vector potential A(f) must be added in the action. In a 
ring of perimeter L pierced by a magnetic flux c/>, we have 
2eA — 9/L, where 9 = 47r0/c/>o is the reduced flux. 



The cooperon has the structure V(r,f;t) = 
V(r,r;t) (e'*)y^ where V(r, f;t) is the cooperon in the 
absence of the electron-electron interaction and (■ • -)c de- 
notes averaging over closed Brownian curves. 

The average over the Gaussian fluctuations of the field 
in l|18l) can be performed thanks to the relation (e'*)v — 
e -3<* >v t erm that appears in the action is 

2e 2 T '* 
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dTW(r(T),r(f)) 



(20) 



where we have introduced the function W, defined in the 
most symmetric way by : 

W(f, f ') = P^r) + P d {r',r') _ p ^ fl) ^ 

For one-dimensional wires of section S we have 
Pd(r,r') — > ^Pd{x, x 1 ), where Pd{x,x') is now the one 
dimensional diffuson. Similar subsitutions holds for W 
and V . The cooperon finally reads 



V(x, x; t) 



x(t)—x 



x(0)=x 
Jo*dr 



Vx{t) 



(22) 



xc 



Mr)' , 2 
4 ° + VTJ T ?/ 2 



W(x(t),x(t)) 



where we have used the definition of the Nyquist time 
(JHJ • The weak localization correction is now given by 



(Ad) = -2 



e 2 D 



dtV{x,x;t) e 



(23) 



which is another way to write eq. J7J) (note that for a 
translation device such as a wire or a ring, V{x,x]t) is 
independent oni). 

Scaling.— If we introduce the dimensionless variables 
u = r/t and y = x/^/Dt, using the expression of Wix, x') 
given below by eq. (|27J) . we get, for a ring or for a finite 
connected wire : 



V(x, x; t) 
-Jo 1 ** 



1 



v(X)=v 



Vy{u) 



Dt Jy(0)=y 

ly(u) 2 +2(^f^\y(u)-v(u)\(l-(^) 1 ' 2 \y(u)-y{u)\) 



(24) 



where u = 1 — u. We obtain the structure assumed in 
the introduction : 



V{x,x;i) = V(x,x;t) x / 



(-,-) 

\TN T D J 



(25) 



where f(x,y) is a dimensionless function. Finally it is 
clear that the integration over time, eq. leads to a 

conductivity of the form 



e- T ILL 
(Act) = -Z x </( 



(26) 
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where g(x,y) is a dimensionless function. 

In the following we will omit the section of the wire. 

How to get rid of time-nonlocality ?— The expres- 
sions H18I22(I were the starting point of ref4 in which the 
correction was computed in the case of an infinite plane 
and an infinite wire. A first difficulty to evaluate the 
path integral (|22|l is the nonlocality in time of the action. 
This problem can be overcome thanks to the translation 
invariance which makes the function W(x, x') a function 
of the difference x — x' only. Such a property is true 
only in few cases. More precisely for the infinite wire 
(AAK) and a finite isolated wire, the function is given 
by W(x,x') — ^\x — x'\. For the connected wire and the 
isolated ring, it reads (see appendix IaT) 



W(x,x') = ~\x- 



1 - 



L 



(27) 



III. EXACT CALCULATION OF THE PATH 
INTEGRAL IN THE ISOLATED RING 

From eqs. Ij22l29|l we can write the weak localization 
in the form : 

2 roo r x(t)=0 

(Act) = -2— / dte~ t/L * / Vx(t) 

n JO Jx(0)=0 



2e^ 

7T 



xe' 



[x = 0\ 



/„' dr \-* r +2iexA{x)-^ r W(x(r),0) 



rV- D * 



(30) 
1 3 = 0) (31) 



where we have rescaled the time to get rid of the diffusion 
constant, and introduced the coupling to the magnetic 
field. Inside the ring the vector potential is given by 
2eA(x) = 8/L, therefore the covariant derivative is D x = 
d/dx-iQ/L. 



The result for L u 



For a translation invariant problem we can follow the 
strategy of AAK : separate the path integral into two 
parts over the time intervals [0, t/2] and [i/2,t] with 

$-Dx(t) - Jdx'J^Vx^J^Vx^r), then 
perform the change of variables R(t) — ^^Hf^ 7 ") anc i 
p ^ = ^i(r)-| 2 (f) ^ the Jacobian is iy Since w(x,x') is 

only function of x — x 1 , the action possesses a potential 
term W(xi(t), X2{f )) = W(p(r),0) function of p only, 
therefore local in time. 

This result is actually due to a general property men- 
tioned in ref. 22 . The path integral (|22|l extends over all 
Brownian paths coming back to their initial value (these 
paths are called Brownian bridges). If we consider a 
Brownian bridge (x(t),0 ^ r ^ t\x(Q) = x(t) = 0) we 
can write the following equality in law : 



x(t) — x(t — t) = x(2t) 



for t e [0,t/2] 



(28) 



As we have mentioned in the introduction, the most 
striking effect of the electron-electron interaction is a 
modification of the dependence of the AAS oscillations 
as a function of L/Ln- For this reason we begin this 
section by emphasizing the results for two limiting cases 
demonstrated below. 

• Small perimeter L <C Ln- In this case : 

{Aan) ^^±lV%( ^V /2 e" 1 " 1 ^^ 372 . (32) 



7T 



L 



• Large perimeter L ^> Ln. This regime is the one stud- 
ied by LMiLi£. We obtain 

(^)**L„m(4)"e- M *W''> (33) 



^ Ln A7(0) 



where - r (V 3 ) 

wneie Ai , (Q) - 3 i/3 r(2 / 3) 



T.372. The exponential 



behaviour coincides exactly with the one obtained by LM 
(see&). However our prefactor differs as will be discussed 
in section Ivl 



(Two random variables distributed according to the 
same probability distribution are said to be "equal in 

law"). Therefore it follows that J * dr V(x(t) — x(f)) (1 =' 
J * dr V(x(t)) for any function V(x). This immediately 
gives 



x(0)=x 

r x(t)=0 



lx(0)=0 



Vx(T) e 



(29) 



B. Computation of the cooperon 

Starting from eq. I|31(l we introduce the rescaled vari- 
able x = x/L. Then the weak localization correction 
rewrites 

e 2 

(Act) = -2— LC(0,0) , (34) 
where the Green's function is solution of 

2 r3 r2~ 



L 2 



c(x,x') 

S(X-X') (35) 
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for x S [0, 1] and periodic boundary conditions. We in- 
troduce the notation 



L 

Ln 



iil id h = | — - 

L, n 



(36) 



We first consider the Cauchy problem. The differential 
equation 



fix) = 



(37) 



is a hypergeometric equation and can be solved by stan- 
dard methods 2 ^ which suggest the following transforma- 
tion f(x) = e~2 s y(s), where the variable is 



e^a 1/4 I x 



(38) 



It is now easy to see that the function y(s) is solution of 
the Hermite equation. A solution of (|37|l is 



f(x)=e-* 8 H„{a), 



(39) 



where H u (s) is the Hermite function (see appendix IB1 . 
The index v reads 



1 

v = 

2 



with oj = 



2^ 



(40) 



Thanks to the symmetry of the differential equation with 
respect to the substitution x *-* 1 — Xi an d since f(x) is 
not invariant under this transformation, another possible 
solution is f(l — x) = e~2 s H y (—s). 

In order to construct the Green's function C(x,x') 
we introduce another solution of eq. satisfying the 
boundary conditions 



/(0) = 1 and /(1) = 0. 
This solution is related to f(x) by 

f(0)f(x) - - X) 



fix) 



/(o) 2 -/(i) 2 



(41) 



(42) 



The derivatives of f(x) at x — an d 1 will be central 
quantities in the following. They can be related to the 
derivatives of f(x) as 



/'(0) = 
/'(I) = 



/(0)/'(0) + /(!)/'(!) 

/(0)2 _ ; (1)2 

/(0)/'(l) + /(l)/'(0) 



/(o) 2 -/(i) 2 

Then the cooperon is given bjsi 2 * 2 ^ : 

e ie(x-x') 



(43) 
(44) 



c(x,x' 



M 



f(x)f(x') + e i0 f(x)f(i-x') 

+e~ w /(l - x ) fix') + /(i - x) /(i - x') 
e^Cx-x') 



/'(I) 



-/(max ( X , x')) /(I " min (X. x')) - (45) 



where M. - 

obtain C(0,0) = i, therefore 



2/'(0) + 2cos0/'(l). At the origin we 
_ j_ 

" M ' 



(Aa) = --L 



1 



7T -/'(0) + /'(!) cos ( 



(46) 



This result shows that in the most general case with ex- 
ponential relaxation (t v ) and electron-electron interac- 
tion (tat), the flux dependence of the weak localization 
correction has still the same structure as AAS, eq. J5J. 
As a consequence the harmonics still decay exponentially 
with n : 



(A(7„) = L 



* V/'(o) 2 -/'(i) 2 ' 



(47) 



where, by analogy with eq. ©, we have introduced an 
effective "perimeter" £ c g, defined by 



cosh- 



eft 



/'(0) 
/'(I) 



(48) 



The equations I|47I48|) are central results of the paper. 
We now analyze several limiting cases, which requires a 
detailed study of /(0), /'(0), /(l) and /'(l). 



C. No electron-electron interaction : Ln = oo. 

We check first that we recover the result © of AAS. In 
terms of the parameters a and b, the limit Ln — ► oo cor- 
responds to a — > with 6 finite, then a; ~ bj (2-y/a) — ► oo. 
With the help of the asymptotic behaviour (|B4f) we 
get f(x) <x e-^x which gives /(x) = 

The derivatives are f(0) = -Vhcothv^ and /'(l) = 
— v6/sinh \/& and the effective perimeter reads £ c ff = 
Vo = L/L v . The harmonics are given by eq. ©. 

D. Small perimeter L <C Ln- 

Instead of expanding the exact solution given above for 
a <C 1, we go back to the differential equation (|37|l and 
construct the solution /(x) by a perturbative approach 
for the small parameter a. This perturbative method is 
explained in the appendix [D] It follows from the expres- 
sions (|D5ID6jl that 

cosh L s = cosh Vb-\ sinh \fb + 0(a 2 ) (49) 

12Vb 

and 

/'(0) 2 - /'(l) 2 = b + JL f C oth ^ - ^) + O(a 2 )(50) 

In the limit L <C the effective perimeter can be 
written 




(51) 
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This combination is not surprising from eq. (|37|l (1/6 is 
the average of %(1 — x) over the interval). The prefactor 
of the harmonic is given by : 



Vfw - far - 4ff 



(52) 



In the opposite limit L 3> L v , the effective perimeter 
can be expanded as 



4ff 



L 



1 L V L 2 



12 L"pj 



and the prefactor is given by 



V/'(o) 2 -/'(i) 




(54) 



For L v — oo, it is clear that eqs. (|51I52|I with eq. 147|) 
give eq. (l3~2*)l . 



r 3/2 
r __kiV_ f 



[Semiclassical i 


[Perturbative \ 

, ul^ 1 y 

! t /g) 


1 % / L \ 3/2 

! 8vl w / ; 






1 _L 




^gS&^_= \ 1 


r"" -Lcp 
i 





function is reduced by an exponential factor e 
this limit u ~ ^ > 1). It follows from that 



Ai'(i^/^) 



f(0) 

/(0) ^ Ai(2^/L») 



(55) 



(this expression is also derived in appendix^lby a differ- 
ent method). 

The relation (|B14)l shows that /'(l) ~ -ie"™/'^). 



(53) Therefore we expect that /'(l)//'(0) - e~™ which leads 
to £ e f{ ~ 7rw. However this dominant term in /'(l) is 
imaginary and does not contribute to /'(l) which is given 
by the next term in the expansion of /'(l). Instead of 
performing a systematic expansion of /'(l), we use the 
semiclassical solution for the cooperon (see appendix Q 
which leads to the behaviour (|C22IC25(1 . As a result the 
effective perimeter is given by the sum of two contribu- 
tions : 



4ff = [-r- 



L 

Ln 



3/2 



where 



r)[x) 



T 3 



arctan 



<L 2 N /Ll) 



(56) 



4.r 



^ (57) 



(see appendix ■ The second term involves the small 
and smooth function 



«(A) = In ( -47ret A3/2 Ai(A)Ai'(A) 



(58) 



k(A) interpolates between k(oo) = at large A and 
k(0) = ln(2/ v y 3) ^ 0.1438 at A = (see figure©. Since 
the first term in i|56|) is much larger than 1 for the regime 
considered in this subsection, /{(L^/L 2 ) can be neglected 



m most cases. 



FIG. 1: (Color online) On this figure we summarize the dif- 
ferent limits for the effective perimeter £ e s given by eq. \4^ - 
For the regime L <C Ln , we have used a perturbative method 
in section [ill D\ The regime L 3> Ln has been studied with 
the semiclassical approximation (instanton approach) in sec- 
tion 1111 M Since the effective perimeter has the structure : 
£cS = (L/ Ljv) 3 ^ 2 r){L? a j L? v ), the line L v — L c oc L 3 J 2 sep- 
arates the regimes of large and small L v . The dashed area 
corresponds to the crossover region where the full expression 
1561) is needed. 



< 



0.15 



0.10 



0.05 



0.00 




-3 -2 -1 

log 10 (A) 



E. 



Large perimeter L 3> Ln. 



In this limit, the function l|39l) . f{x) = 
e -i±VS( x -i/2) 2 i/_i +iM ( e W4 a i/4(^ _ 1/2)), presents for 

X — > a behaviour given by l|B10|l . For x ^ 1 the 



FIG. 2: The function k(A) of eq. JKfy. 



Finally the weak localization correction reads 



7T 



^'(L%/Ll) 



(59) 
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We remark that the harmonic coincides, as it should 
in the limit L ^> Ln, with the result of AAK, eq. 
for the infinite wire. This provides another check of the 
exact solution, and in particular of its prefactor. 

We now discuss the various limiting cases obtained 
by varying L v . First we remark that the prefactor of 
eq. lj5"9"|) has the form L^ g\{L 2 N j L 2 ^), where 51(2;) is 
a dimensionless function, whereas the effective perime- 
ter has the form £ eS = (L/L N ) 3 / 2 n(L 2 c /L 2 v ) with L c = 

L 7v \J Ln J ' L <C Ln (if we neglect the smooth contribu- 
tion). Therefore we have to distinguish three regimes : 

• Negligible exponential dephasing : L v 3> Lat. 
Eqs. 1561591) give eq. H33[) . This is the only regime in 
which the contribution k(A) in eq. (|56|l plays a role. 



• L c — Ln\/ Ln / L <C L v <C Ln- 

fies as 



The prefactor simpli- 



(At7„ 



T e -|n|4ff 



and the effective perimeter can be expanded as 



4ff 



TT / L 

8 [J^ 



3/2 tt L^L^ Ail 

+ 2 l% : 



3/2 



(60) 



(61) 



The two first terms correspond to itlo. The effective 
perimeter is dominated by the first term. However the 
function 77(6/0) appears in the argument of an expo- 
nential, in eq. (|59|l . multiplied by a large parameter as 
4ff = y/ar)(b/a). Therefore it is not clear a priori when 
the terms of the expansion of r](x) are negligible. 

• Dominant exponential dephasing : L v <C L c . In this 
case the expansion of the effective perimeter reads 



4ff 



L 

Lin 



1 L V L 2 
12 L% 



(62) 



which coincides with the expansion of the result £ c ff — 
\Jb + a/6 ~ \fb + YiVS obtained by a perturbative ex- 
pansion in a. If L p <C L^{Lm / ' L) 2 , we recover the result 
of AAS ©■ 

We have summarized the different limits for the effec- 
tive perimeter on figure ^ 



IV. RELAXATION OF PHASE COHERENCE 

In this section we interpret the results of the previ- 
ous section in a time representation and give a rigorous 
presentation of the heuristic discussion of the introduc- 
tion. The results of this section may be useful to consider 
more complicated situations than an isolated ring, when 
the path integral cannot be computed exactly, like the 
connected ring studied in the next section. Let us con- 
sider the n-th harmonic V n (x, x' '; t) of the cooperon. The 



Fourier transform over the magnetic flux of the path in- 
tegrals, eqs. (|18I22|) in which the coupling to the external 
magnetic field has been re-introduced, selects the paths 
with a winding number equal to n. We can write : 



P n (x,x;t) 



fx{t)—x 



lx(0)=x 



2tt 



—V{x,x;t) e 



And 



(63) 



where A/"[x(r)] is the winding number of the trajectory 
around the ring. For a closed trajectory (x(0) — x{t)) 
we have A[[x{t)] = j- / *drx(r). (In this section we set 
D = 1). Let us introduce the probability V n (x,x';t) for 
a Brownian curve to go from x' to a; in a time t encircling 
n times the flux : 

V n (x, x'-t)= ( { ) Vx(t) e-Ho d-M 2 S nM[x(T)] (64) 

Jx(0)=x' 



For an isolated ring this probability simply reads 

V n (x, x';t) = -= e- <X " X «" L) " . (65) 



Then we can rewrite the harmonics of the conductivity 



T^n(x 7 X,t^ — 7^7}(x,X,£) {q 



v,c n 



= V n (x,x;t) (e-^v) 

with ($ 2 )y given by (|2UJ| : 
1 2 '* 



$2N 



3/2 , 



&tW{x{t),x{t)) 



The average 



■ c n 



1 

V n {x,x;t) 

fX{t)—X 

■ i Vx(t) ■ ■ ■ 

c(Q)=x 



(66) 
(67) 

(68) 

(69) 



l/o dTi ( T ) 2 s. 



n,N[x(T)] 



is performed over all closed Brownian trajectories with 
winding n. The probability V n in the denominator en- 
sures the normalization (l)c„ = 1< Therefore the func- 
tion (e 1<s> )v,CnJ related to the n-th harmonic of the AAS 
oscillations, characterizes the relaxation of phase coher- 
ence due to the electron-electron interaction for trajecto- 
ries with winding number n. We now analyze this quan- 
tity. 



A. Diffusion of the phase 

An indication on the nature of the phase relaxation, 
characterized by the function (e _ 5(* ' )v )c n , can be ob- 
tained by studying the diffusion of the phase, that is the 
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much simpler quantity (<& 2 )v,C„- If we consider r < f/2, 
the average (W)c n is given by 

V n (x,x;t)(W(x(T),x(f))) c =y [ dxdx' W(x, x') 

L Jo 

+00 

(x>,x;t~2r).(70) 

m— — 00 

The final result is symmetric with respect to 2r «-> £ — 
2r. The double integration can be reduced to a simple 
integration thanks to the relation L dxdx' f(x — x') — 

Jq du (1 — u)[f(u)+f(— u)]. Then the integral is unfolded 
to extend over R. We obtain : 

/+00 
dxn(x) 
-00 

( a (x-^-L) 2 \ 

e _ sF e 4 c- 2t ) \ 
v 87tt ^/47r(i - 2t) / 

where Q,(x) is the function Q(x) = 2W(x,0)(l - f ) = 
x(l— ^) 2 for x G [0, L] and periodised on R. To go further 
we distinguish two cases depending on the relative order 
of magnitude of the time t and the Thouless time 773. 



The integration over time leads to average the function 
n : f* dt (W) Cn = t / L ^ It immediately follows 

that : 

N c 

Ergodic regime (t 3> tb j 

The cases n — and n ^ can be treated on the 
same footing. The Gaussian function in eq. (|71|l is very 
broad compared to L and we can replace fl(x) by its 
average value Q(x) — > j, Jo dxil(x) = It follows 
that (W) Cn ^ ^, then : 

>Vc„4^ Vn. (76) 

To summarize we see that, for the harmonic n — 0, the 
diffusion of the phase crosses over from a i 3 / 2 behaviour 
to a linear t behaviour, whereas for n 7^ it behaves 
always linearly. This difference shows up in the function 
(e 1$ )y i c„ leading either to a non exponential or to an 
exponential phase coherence relaxation. 



Diffusive regime (t <C To ) 

In this case we have to separate the cases n = and 
n ^ 0. 

• Harmonic n — 0. The integral is dominated by the 
neighbourhood of x ~ since the Gaussian function is 
very narrow compared to L. We can replace the function 
Q(x) by its behaviour near the origin : £l(x) — > 9{x)x, 
where 8(x) is the Heaviside function. We obtain 



(W(x(r),x(f))) Co 



2r(t - 2t) 

7rf 



(72) 



This result should be symmetrised for r > t/2. Integra- 
tion over time r gives 



1 



4> 2 



'V,C 



3/2 



(73) 



In this regime the geometry plays no role and we recover 
the result obtained for an infinite wir o 10 i 13 . This result 
is related to the AAK behaviour for ijy 3> L v . 

• Harmonic n ^ 0. For times t < tdi the integral 
(|71|l can be estimated by the steepest descent method : 
it is dominated by the neighbourhood of the point mini- 

- \2 

that is x = nL2r/t. Therefore we 



mizmg =- 

or 

obtain : 



(x- 



4(*-2t) 



(W(x(t),x(t))) c 



[tt{nL2T/t) + Q(-nL2r/t)] (.74) 



B. The function (e i *) v ,c„ 

The calculation of (Q 1 ^)v,c n is a more difficult task. It 
can be obtained by different strategies. 

(A) Small phase approximation.- At short times, the 
phase $ is small, therefore we can linearize the expo- 
nential so that : (e~i^ v ) Cn - e"^* 2>v '- c ™_Then we 
can use the results given above in subsection IIV Al 

(B) Inverse Laplace transform.- The weak local- 
ization correction to the conductivity (Acr n ) ~ 
/ °° dt e - */ 7 ''' V n (x, x; t) (e 1 *)y,c„ has been derived for ar- 
bitrary t v — L 2 . Physically, the parameter t v takes into 
account other dephasing mechanisms responsible for an 
exponential relaxation of phase coherence. From a tech- 
nical point of view the parameter t v allows to probe time 
scale t ~ t v in the path integral and we can, in principle, 
compute the inverse Laplace tranform of (Acr n ). 

(C) Large phase for n 7^ 0.- If none of the previous meth- 
ods can be used (the method (B) because it is too diffi- 
cult, and the method (A) because it is not the range of 
interest), we can use the following remark : for t v — 00 
we see from (|32|) that 



(Acr„) 



, 1 -mz£ 
dt —=e 4t 

y/i 



(77) 



We expect that the behaviour at large time involves 
the tail of (e 1 *)y.c„ which we can assume to behave as 
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(e i *)v,c„ oc £^e _Tt . The integral of the l.h.s can be 
estimated by the steepest descent method 



, 1 
d£ —=e 

y/i 



-rt a 



2- 



a(a + l)T 



2+2 „-T(a+l)t; 



(78) 



where t* = [^f^) 2 ]^- The coefficient T and the 
exponents a and /i are obtained by comparison of the 
dependence of this result with L and n with the known 
behaviour for (Act,,). 



C. Small perimeter L <C Ln 

We now analyze (e 1 *)y j c„ in the small perimeter limit. 

a. Short times 

In the short time limit, the linearization of the ex- 
ponential is valid (method (A)). Therefore we can use 
expressions (|73I75I76|) . These expressions give a pre- 
cise definition of the "short time" regime, which extends 
until ( ( E >2 )y i c 7l ~ 1 that is £ ~ r c . The time scale 
r c , given by eq. l(Tfi]l. is associated to the length scale 
L c = Ln^/ITnJL introduced in subsection IIII El fsee^. 



c. Summary 

From all these results we can conclude that for har- 
monic n = the relaxation is non exponential at very 
short times and eventually becomes exponential for time 
larger than the Thouless time : 

(e^y.Co * exp-^f— Y for£«r Z5 (82) 
4 V T N J 



exp - 



1_£_ 

6 Tc 



for £>ri5. (83) 



This difference comes from the time evolution of (W) c : 
when the ring has not been explored (£ -C Tjj) it scales 
like (W) Co ~ \fi, while it becomes time independent for 
ergodic regime (£ ^> td). 

On the other hand the phase coherence relaxation is 
always exponential for harmonic n =/= : 



exp - 



1 t 

6 Tc 



V£ 



(84) 



This is due to the fact that the trajectories with fi- 
nite winding necessarily explore the ring which leads to 



( W )c n 
tion. 



£ , for all times, as explained in the introduc- 



D. Large perimeter L Li 



b. Long times 

In this case we consider the harmonic n = and n =/= 
on the same footing. The regime £ S> tjj corresponds 
to L <C hu>. With L <C Ln this leads to the "pertur- 
bative" regime a,b <ti 1 for which we have found the 
expressions I|51I52[) : 



/>oo 

/ d£e~* /r " T n (x,x;t) 
Jo 




The inverse Laplace transform can be computed exactly 
in this case. It gives 

V n (x,x;t) = —j=e ~e . (80) 
2v 7r£ 

We immediately obtain 

( e ^)v.c = CX P ~\— for f » T o V,i - ( 81 ) 

This is the same result as for i<r c . Since tjj <C r c , it 
turns out that the result obtained from the linearization 
of the exponential is valid for all times. 



a. Short times 

In this regime the analysis provided for small perime- 
ter using the results of section llV Al remains valid. The 
condition of validity of the results slightly changes since 
the times are now in the following order : r c <C tjv <C td- 
For harmonic n = 0, eq. (|82|l holds for £ <C t^. For the 
harmonics n ^ 0, eq. I|84|l holds for £ <C t c . 

b. Long times 

• Harmonic n — 0. For this harmonic, the mag- 
netoconductivity is given by AAK since for Ljv < I 
the boundary conditions are not important. The inverse 
Laplace transform of eq. @ has been computed exactly 

inref3i^: (e^) v , Co = ^En=i 
where «„'s are zeros of Ai'(z). 

• Harmonic n^O. In this case we can only use the 
method (C). We compare to the integral GHJ : the 
dependence in n of the exponential gives the exponent 
a = 1. Then its L-dependence gives T = ^-h-- The 
analysis of the prefactor shows that \x = 0. Therefore : 

(e^c^exp-jj* (85) 
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c. Summary 
For the harmonic OiLii : 

(e^) v . Co =s exp-^f— ) 7 fart*™ (86) 
7r£ 1 i 

— 1 — : exp- ui — lor tat < t (87) 

T N \Ul\ T N 

(the first zero of Ai'(z) is \m\ ~ 1.019). 
For harmonic n^O, we have seen above that : 

(e i<E> )v, c „ ^ exp-i— fort<T c (88) 

IT 2 t 

~ exp for r c < t . (89) 

64 r c 



E. From exponential phase coherence relaxation to 
non exponential size dependent harmonics 



SemicTassTcal 



\nt 1 ii t 
x-,IN eX P-' Mll % 




Perturbative 



exp- 



6i c 



fktt \ 3/2 



FIG. 3: (Color online) Relaxation of phase coherence for tra- 
jectories that do not wind around the flux : (& }v,c - The 
result for the infinite wire (AAK) is recovered when trajecto- 
ries cannot explore the whole ring, that is either for t -C td 
or for t <C tjv . 



On figures and 0] we summarize the results obtained 
for the function {e l ^)v,C n - 

The behaviour ($ 2 )y oc i 3 / 2 was first mentioned in 
ref^, where it was conjectured that it may lead to inter- 
esting effects in a ring. However, when the effect of wind- 
ing is properly taken into account, it turns out that the 
interesting effects in the ring come from an exponential 
relaxation, i.e. )v,c n <x t. In order to emphasize this 
point, let us summarize the relationship between time 
dependence of the phase relaxation and the decay of the 
harmonics. For n ^ 0, the function (B )v,C n is always 
exponential, exp— /3t/r c , with f3 = 1/6 or (3 — 7r 2 /64, 
depending on the time regime (see figure 0}. The weak 
localization is given by the time integrated probability to 
turn n times around the ring weighted by the exponential 
damping : 




(Aoy, 



1 



y/t 



(90) 



We recover the non exponential size dependence of ( Acr n ) , 
eqs. IM . 



FIG. 4: (Color online) Relaxation of phase coherence for tra- 
jectories with a finite winding number : (6 1 )v;e„) unth 



(Aa n ) 



exp - 



L \ ^ -\n\L 3 / 2 T^ 2 



J N 



3/2 



(91) 



consequence of an exponential relaxation of phase coher- 



V. THE EFFECT OF CONNECTING WIRES 

Up to now we have considered an isolated ring. This 
was an important assumption in order to calculate the 
path integral. However, in a transport experiment the 
ring is necessarily connected to wires through which the 
current is injected. This has two important consequences 
that we now discuss. 



A. Classical nonlocality and quantum nonlocality 



• Classical nonlocality.- The classical conductance of a 
wire of section S and length L is given by the Ohm's law 
G cl = ctqS/L where a is the Drude conductivity. This 
result can be rewritten for the dimensionless conductance 
as g^ irc = 2§/x = a d N c £ e /L, where iV c is the number of 
channels, l e the elastic mean free path and a<j a numerical 



constant depending on the dimension {a.\ = 2, ct2 = 7r/2 
and «3 = 4/3). The quantum correction to the classical 
result is given by 



2 f L 

(Ag wilc ) = -—J dxP c (x,x), 



(92) 
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where we have introduced the notation P c (x,x') = 
f™ dte- t / T vV(x, x';t). 

For a multiterminal network with arbitrary topology, 
the classical transport is described by a conductance ma- 
trix that can be obtained by classical Kirchhoff laws. 
This classical conductance matrix is a nonlocal object 
since each matrix element depends on the whole network 
and the way it is connected to external contacts. On such 
a network, because of the absence of translation invari- 
ance, we have shown in ref^i how the cooperon must be 
properly weighted when integrated over the wires of the 
network in order to get the weak localization correction 
to the conductance matrix elements. Eq. (|92|) generalizes 
as a sum of contributions of the different wires : 



2 y-^dC 



dU 



dxP c (x,x) , 



(93) 



where C is the equivalent length obtained from Kirch- 
hoff laws : g cl — adN c £ e /£. The weight of the wire i 
is the derivative of the equivalent length with respect to 
the length of the wire k. The existence of these weights 
can lead to unexpected results, like a change in sign of 
the weak localization correction for multiterminal geome- 



When we consider the ring of figure El the equivalent 
length is given by C = l a + l c \\d + h where l^ 1 
It follows that 



- r 1 +r 1 



<A 5 ) = - 



('o + l c \\d + h) 2 
I 2 , f I 2 



a (lc + ld) 2 J c (lc+ld) 2 J d 



(94) 

da; P c {x, x) . 



• Quantum nonlocality of the cooperon.- The cooperon 
is a nonlocal object that depends on the whole network. 
P c {x, x) is a sum of contributions of diffusive loops that 
explore the network over distances of order of the phase 
coherence length. We have shown recently in refsiSLSS 
that the presence of the connecting wires can strongly 
affect the behaviour of the harmonics of the AAS oscilla- 
tions. We can distinguish two regimes for long connecting 
wires : (i) in the limit L v <C L the AAS harmonics are 
exponential (Ag n ) oc exp — \n\L/L v . (ii) However in the 
limit L <C Lip, the behaviour of the harmonics becomes 
(Ag„) oc exp — \n\ y j / 2L/ L v . This different behaviour was 
analyzed in detail and shown to originate from the fact 
that the Brownian trajectories can explore the connecting 
wires over distances larger than the perimeter L. The ef- 
fective perimeter of a Brownian trajectory encircling the 
ring is £ e ff — y/2L/L v » L/L v . 

We see that the simple exponential decay of the har- 
monics with the perimeter, eq. JSJl, can be modified for 
two reasons : either the presence of connecting wires, or 
the effect of the electron-electron interaction. One acts 
on the nature of the diffusion around the ring, the other 
acts on the nature of the dephasing. In this section we 
propose to combine these two effects. 



The presence of connecting wires modifies the diffuson, 
and therefore the function W(x, x'). The main difficulty 
to compute the path integral giving the cooperon is that 
W(x,x') can not be written as a single function of x — x' 
and it is not possible to make the path integral local in 
time. 



a 




FIG. 5: (Color online) A mesoscopic ring connected at two 
reservoirs (represented by the wavy lines). 



B. Large perimeter L 3> Ljy 

The connecting wires are not expected to have a strik- 
ing effect in this regime. Since the cooperon vanishes ex- 
ponentially over distances larger than Ln (or L v ), the 
integration of the cooperon over the connecting wires 
can be neglected when studying the harmonics. On the 
other hand the diffuson is affected by the presence of 
connecting wires, and the function l|27|l has to be re- 
placed by l| A4I A6|l . This function has a similar struc- 
ture to l|27|l ; moreover the two functions are equal in 
the limit of long connecting wires as discussed in ap- 
pendix The function was given in refill for a sym- 
metric ring (l a = lb and l c = Id). In this case, when 
both x and x' are in the same arm of the ring it readsll 
W(x,x') = \\x-x'\{l-^\x-x'\), where 7 = l c \\ d /C. 
This explains why, in the regime L 3> Ljv> the pres- 
ence of connecting wires has almost no effect. It es- 
sentially modifies the numerical prefactor in the expo- 
nential : (A<7„) cx exp — n^-|(-^) 3 / 2 . The coefficient 
Cj/Cq interpolates smoothly between 1 (for 7 = 0) and 
1/V2 (6*7=1)11. 

Prefactor of the harmonics. - Before going to the limit of 
small perimeter we consider into more details the prefac- 
tor of the conductance. The exact result l|59|l was derived 
for an isolated ring and it is not clear how it is related 
to the conductance through a ring connected to contacts 
by arms. In this latter case the conductance is given by 
eq. B94JI . In the limit Ln <C la, h -C l c ,ld, let us approx- 
imate the cooperon by 



P c (x,x) 



P c \ 00 wire for x £ a,b (95) 

-P C | isolated ring for X € C, d , (96) 



where P c |oo wirc and P c | isolated ring are the cooperon for 
an infinite wire and the isolated ring, respectively. This 
means that we neglect the fact that the cooperon can 
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leak in the wires a and b, when its coordinate is inside 
the ring. The harmonic n — 0, given by 



L N M(L%/Ll) 



(97) 



has the form (Ago) = (w '/ 'e 2 )(A<7o) / The harmonics 
ti^O read 



(Aft, 



i C ||d 



£ 2 A\{L 2 N /L 2 ) 



(98) 



where the effective perimeter is given by eq. i|56[l. 

One may question the validity of the hypothesis 
I|95I96[I . To answer this question, let us consider the limit 
Ln/Lu, = oo of eqs. I|97I98|) and compare with the exact 
result obtained in ref. 12 . In eqs. I|97I98|I . the ratio of Airy 
functions is replaced by —L v /Ln and £ e ff = LfL v . On 



the other hand the exact result gives (Ago) 
and, for the harmonics n ^ 0, 



<A<7„) 



c\\d n tp 



2|n| 



-\n\L/L v 



-L v /C 



(99) 



and eq. (JSHJ for 

missed the factor 



The comparison between eq. (|99|) 
Ln/L v — oo shows that we only 
(2/3) 2 '™'. This factor is related to the probability for the 
diffusive trajectory to remain inside the ring, when arriv- 
ing at the vertex (see refsii 2 ^ 2 ^* 3 ^). Note that, in the limit 
of short connecting wires l a <C L v , the factor (2/3) 2 '™' in 
eq. © is replaced by [1 + L 2 /(Z Q L)](2Z a AM N , which 
describes the probability that the diffusive particule is 
not absorbed by the nearby reservoir when arriving at 
the vertex. In the limit l a /L v — > 0, the reservoirs break 
phase coherence at the vertices and the harmonics vanish. 

Finally we remark that the prefactor of the harmonics 
differs from the one obtained by LMi£. With our nota- 
tions their result reads : (Ag„)LM ~ L 9 ^ 4 e"'™'^ 17 " ■* , 
whereas our prefactor is linear in Ln (for L v — oo), as 
for (A<7 n ). We stress that the difference between LM's 
and our result is not due to the presence of connecting 
wires. The linear dependence in Ln of the prefactors 
of eqs. H59I98|I comes from the path integral. In LM's 
paper, the L 9 ^ 4 comes from a wrong estimation of the 
prefactor of the path integral. In appendix [U] we have 
shown how the correct prefactor can be extracted within 
the instanton approach followed by LM. 



C. Small perimeter L <C Ljv 

In this limit, the Brownian trajectories contributing to 
the path integral (|18I22|I are related to times t ^> td . It is 
known that for such time scales the arms have a striking 
effect since the diffusive trajectories spend most of the 
time in the long connecting wires (see ref^i and section 



5.5 of refi2£). This affects both the winding properties 
around the ring and the nature of the dephasing. 

We expect that the relaxation of the phase coherence 
mainly occurs inside the arms, i.e. the largest contri- 
bution to (3> 2 )y oc J Q dTW(x(T),x(f)) correspond to x 
and x' in the arms. In this case the function W(x, x') is 
given by eqs. I|A7IA8(I . If we consider the long arm limit 
la,h 3> Ln ^ hJd, we can take the limit l a ,h — * 00 
in (|A7IA8p and we have : W(x I. We re- 

cover the same function as for the infinite wire. In this 
limit the length over which the trajectories extend in the 
wires is not limited by the size of the system but by the 
time. Therefore we expect that the phase coherence re- 
laxation is non exponential and the function (e 1 *)y i c„ is 
similar to the one obtained for the infinite wire (or for 
the harmonic n = for the large ring, eqs. 1)861871) ) : 
(e i *)y, c „ ~ (e i *)y, c |wire = /wto(*/rv). This funtion de- 
creases monotoneously over the scale tn- 

Since the cooperon P c (x,x) is expected to decay over 
a scale Ln in the arms, the integration over x in eq. 1)94)) 
leads to (integration is dominated by the contributions 
of the arms) : 



dtV n {x,x;t){Q^) Vfin (100) 



where x is now a point inside the ring (in the regime 
t » td oc L 2 , the probability V n {x,x;t) is al- 
most independent on x, provided it remains at a dis- 
tance smaller than \ft from the ring). The probabil- 
ity to wind n times around the connected ring is 21 

V n (x,x;t) ~ jj^rrVte = *0) where V(0 ^ 

^(C/4) 1 / 3 e- 3 («/ 4 ) 4/3 for£ > 1 [ip(0) is finite 21 ]. There- 



fore, as a function of time, V n (x,x;t) increases until 
t ~ ti 4 td and then presents a smooth tail i~ 3 / 4 . In or- 
der to evaluate eq. 1)100)1 . we have now to consider two 
limits : 

• For n 2 L <C Ln, the integral 1)100)1 is dominated by the 
tail of the probability V n (x,x;t), therefore 



(Aft 



iff IV* 

(L+k) 2 



(101) 



y/ tj T N exp 

method gives 



For n 2 L 3> Ln, only the tail of /whe^/ 7 ^) ~ 
u\\t/TN is important. The steepest descent 



7-3/2.1/2 / 21- \ 7/12 

^--tt^It^) e " K| * ( 102 ) 

{I a + hV \ L N J 



whith k 



V2\ 



1/4 



1.421. We recover a depen- 



dence of the harmonics reminiscent to the one obtained 
for exponential relaxation of phase coherence 21 - (see dis- 
cussion of paragraph V.A). Finally, we remark that the 
expected decay of the harmonics in temperature is : 
(Ag n ) ~ _T- 11 / 36 e-™ Ll/2Tl/6 , for sufficiently large n. 
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This prediction should be tested experimentally on a 
chain of rings separated by sufficiently long wires, com- 
pared to the phase coherence length (several rings are 
required in order to perform a disorder average). 





Exponential relaxation 




exp —nL/Ltp 
exp-n(2L/L v ) 1/2 




Electron-electron interaction 


L > L N 
L < Liv 


exp-fn(L/Liv) 3/2 
exp — Kn{L/ Ln) 1 ^ 2 



TABLE I: The harmonics {Ag n ) of the conductance through a 
ring of perimeter L connected to long arms. We compare the 
results for exponential phase coherence relaxation, described 
by Ltp, and the one obtained for electron- electron interaction 
(List), k is the dimensionless constant given in the text. 



VI. CONCLUSION 

We have considered the effect of the electron-electron 
interaction on the weak localization correction for a dif- 
fusive ring. We have calculated exactly the path inte- 
gral giving the weak localization correction for the iso- 
lated ring in the presence of electron-electron interac- 
tion (characterized by the Nyquist length Ln) and of 
other dephasing mechanisms described by an exponen- 
tial phase coherence relaxation (characterized by L v ). 
The harmonics of the conductivity are always of the form 
(Aer n ) cx exp —|n|^ c ff , where £ c g accounts for both kinds 
of relaxation, combined in a nontrivial way. The effective 
perimeter can always be written as : 

4ff = (j^)' 2 f{Ll/L%), (103) 

where L c = L 3 J 2 /L 1 / 2 . For large perimeter L 3> Ln, the 
dimensionless function is f(x) — r](x) [eq. (|C26f) ] . For 
small perimeter L <C L N it is given by f(x) = y/l/6 + x. 
All limiting behaviours of £ e ff have been studied in sec- 
tions EDS and EM 

In order to interpret these results, we have studied 
the function (e 1 *)^^ characterizing the phase coherence 
relaxation for trajectories with winding n (involved in 
the n-th harmonic of the AAS oscillations). We have 
shown that, whereas the phase relaxation crosses over 
from a non exponential behaviour to an exponential be- 
haviour for the harmonics n = Q, it is always exponen- 
tial for n =/= 0. The time r c characterizing the expo- 
nential relaxation (e 1 *)v r ,c„ ~ exp— t/r c is given by : 
t c = <ToS/(e 2 TL) = t n ^t d 1 ^ 2 . This exponential relax- 
ation is at the origin of the non exponential decay of the 
harmonics with the size and the new temperature depen- 
dence^ predicted by LM : (Ag n ) ~ ex.p—n(L/L N ) 3 ^ 2 ~ 
cxp-nT^L 3 / 2 . 



In the light of these results it seems difficult to 
interpret the experiment recently performed on large 
GaAs/GaAlAs square networks where a behaviour 
exp— nT 1 / 3 was observedSS^ 3 .. However the experiment 
was performed in a temperature range where Ln ~ L. 
In this regime the diffusive trajectories start to explore 
the network surrounding the loop, which modifies the 
behaviour of the harmonics as a function of L/Ln and 
therefore the temperature dependence. However we have 
not been able to extend the theory to the square network. 

In the last part of our article we have studied the 
effect of the wires connecting the ring to reservoirs. 
Whereas the AAS harmonics are weakly affected by the 
connecting wires in the limit L Ln, we have shown 
that a strong modification is expected in the opposite 
limit L <C Ln , where we have predicted a behaviour 
(Ag n ) ~ T- 11 / 36 exp-nL 1 / 2 T 1 / 6 . The appropriate ex- 
perimental setup to test this result is a chain of rings 
separated by wires whose length remain larger than Ln- 
This situation is particularly interesting because the flux 
sensitivity is due to the motion inside the ring whereas, 
in this case, the dephasing occurs mostly in the arms. 

An interesting effect has been observed recently in the 
study of four terminal measurements of AB oscillations 
in a ballistic ring. It has been shown experimentally 34 
that the dephasing rate depends on the configuration of 
the voltage probes and current probes. It was suggested 
that a measurement with current probes on both sides 
of the ring favors charge fluctuations inside the ring and 
leads to a high dephasing rate, whereas a nonlocal mea- 
surement with current probes at one side and voltage 
probes at the other side (i.e. no current flows through 
the ring on average) diminishes charge fluctuations in the 
ring and therefore leads to smaller dephasing rate. This 
effect has been described theoretically in refw An in- 
teresting question is whether a similar effect might occur 
in a diffusive ring. 
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APPENDIX A: THE FUNCTION W{x,x') 

1. Isolated ring 

We give here the solution of the equation (7 — 
D 2 )P(x 1 x 1 ) = 5(x — x') on a ring pierced by a flux. 
D-r = d/dx — \0/L is the covariant derivative. We 
introduce the variable x = X /L € [0,1]- The solu- 
tion of the equation (b — D 2 ,)C(x,x') = <Hx — x') with 
b = -fL 2 and D x = d/dx — id is given by l|45|l for 
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f(x) — Slnh , , V X ^ . The two Green's functions are re- 

J sinhVfc 

lated by P(x,x') = L C(x/L, x'/L). 

& 6 ^ 0.- Let us consider the limit 7 = 0. We 



• 7 = 

have : 



c(x,x') 



J0(x-X') 



2(1 



(Al) 



1 + lx-x'l (e- [9si ^^ - l) 



• 7 = & = 0.- In the limit of vanishing 7 and 9, 
the diffusion equation possesses a zero mode, therefore 
the Green's function C(x,x') presents a diverging con- 
tribution: C( X ,x') = W ~ \\x~x'\ (1 - Ix - x!\) • This 
diverging contribution disappears when considering the 
function 

W(x,x') = P{X ' X) \ P{X '' XI) -P{x,x') (A2) 



1, nf, \ x ~ x'\ 

= ^\x-x\\l- — — 



(A3) 



The existence of a zero mode is an artefact coming 
from the fact that the system is isolated. In a more re- 
alistic situation (when the ring is connected to reservoirs 
through wires, for example) the Laplace operator does 
not possess a zero mode2&. Physically, the zero mode 
does not contribute to the function W(x,x'), i.e. to the 
dephasing, since it corresponds to uniform fluctuations 
of the electric potential that do not contribute to the 
phase $, given by eq. I|19fl . 



2. Connected ring 

We now construct the symmetric function W(x,x f ), 
defined in eqs. (|21IA2II . when the ring is connected to 



reservoirs by two wires (figure[5j). In this case W(x, x') is 
fully characterized by a set of components corresponding 
to the coordinates in the different wires. A system of 
coordinate must be specified : we first give an orientation 
to the wires of the network, shown on figure (we call 
"arc" an oriented wire). Then the coordinate along an 
arc i belongs to the interval [0, Zj, where U is the length 
of the arc. Below we construct W(x,x') when x and x' 
are both in the ring or both in the arms. 



a 




b 



FIG. 6: Orientation of the arcs of the ring. 
a. Inside the ring 
When both coordinates are in the arc c, we obtain : 
W„{x, x') = \\x- A (l - + 'rif ) ^ 

* \ la + 'c||d + '6 'c + Id J 

with x,x' G [0, l c ]. For l c = Id we recover the expression 
given in refii£. We also need the case when see and 



W c , d (x,x') 



1 (la + h)lc\\d 

2 la + l c \\d + h 
1 " 

"2 



X 1 



X X 
lc Id 

+ x'(l-- 

Ir.J V Id 



,(A5) 



with x £ [0, l c ], x' £ [0, Id] and with the orientation of 
figure El It is more convenient to consider a unique way 
to measure both coordinates. Therefore we shift x' by l c 
in HA5j> : 



W c ^ d (x,x' - l c ) = - 



rt 



la ~t~ lc\\d 
la + ld + h 



la — l c \\d + lb 

X - 

la + Zc||d + lb 

la + lc + lb 
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(lc+ld)(la+l c \\d + lb) 
12 l a + lb ^XX' 



la + l c \\d + h lc + Id la + l c \\d + h lc + Id la + l c \\d + lb l C + Id 



(A6) 



with x £ [0 7 l c ] and x' £ [l c ,lc + Id]- It is now clear 
that in the limit of long connecting wires, l a , lb 3> lc, Id, 
eqs. (IA4IA6II lead to the same result as in the isolated 
ring, eq. I|A3(1 . We see that, inside the ring, W(x,x') is 
not everywhere a function of x — x' only, apart in the 
limit I a , lb — * 00. 



b. Inside the arms 

When both coordinates are in the same arm, we have : 

lV „ ( ,, I -) = i|,-xl(l- i -^L) (A7) 

When the two coordinates belong to different arcs we 
prefer to shift the origin of the coordinate x' , as we did 
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inside the ring. If the shift is chosen to be l a + l c \\di we 
obtain the simple expression : 

W a , b (x, x'-l a - l c]ld ) = Ux' -x)(l- - (Vs) 

(in this expression x' = l a + i c iu corresponds to the 
begining of the arc b). When 1^ = (no ring) we 
obtain from l|A7IA8|) the result for a connected wire 

W(x,x') = — x'\(l — j—n ■ ), which is similar to the 
one for the isolated ring eq. as mentioned above. 

It is remarkable that, in the presence of the ring, there 
exists a choice of coordinates for which W(x, x') in the 
arms has precisely the same structure as in the absence 
of the ring. In the limit of an infinite wire, l a ,h — ► oo, 
we recover from (|A7IA8|) the result of the infinite wire, 
W(x,x') ~ \\x-x'\. 

APPENDIX B: HERMITE FUNCTION 



tion in i|B2|) . One obtains 

2 v / r(i-iw)y >A 

where the phase reads 

(p(x,x) = x 2 + 2Ax + win a; . (B6) 
We introduced the notation 

^ = « 1/4 (x-0- (B7) 

We are interested in the limit a — > oo with L^v / finite 
(or zero), therefore it is convenient to write : 

where A = (Ln/L v ) 2 . 



Consider the Hermite equation 23 : 

y"(z)-2zy'{z) + 2vy(z) = 0. 



(Bl) 



Two linearly independent solutions are the Hermite func- 
tion H v {z) and H v (—z). An integral representation is 



H v {z) = 



1 



dt 



n-v)j f +l 

from which we get the series representation 



H v {z) = 



1 



2T{-v] 



X>1)T 



n — u\ (2z) 



(B2) 



(B3) 



We now study several limiting behaviours of the Her- 
mite function H_ i +ia) (e i7r ' /4 a 1 / 4 (x -1/2)), where u = 

1 I b _r_ \/a\ 

2 ^ yfa i )■ 



The limit a 







In this case u> ~ Jyj' therefore we study the limit 
lu — > oo when the argument of the Hermite function 
reads z = xj \J~\lo with x finite. Using the expression 
2 g-mf ; valid for w 3> n, and the 



if) 



series representation (|B3j) . we get 



H i 



3 -%/2cc 



— VjJ 7 UJ — >oo 



(B4) 



(i) The case x 







The function ip(x, x) is a monotonous function of the 
variable x G M + , however its second derivative vanishes 
at x = ±yJuj/2. For x ~~ * the first derivative at this 
point becomes very small in the limit a — ► oo, therefore 
we expect that the neighbourhood of yJu>/2 brings the 
dominant contribution to the integral. The expansion 
of the phase in the neighbourhood of y/cJ~/2 reads (for 
t -► 0) 

ip(y^]2 + a x ' x %x) = <p(yfc/2,x) (B9) 
+ (a^x + h)t+ l -t i - + ■■■ + O(a^V) . 

Inserting this expression into the integral representa- 
tion (|B5() . we obtain the Airy function 37 Ai(x) — 

i f+™dte i( -£+ xt \ namely : 



H_ 



oc e 

a — *oo 



Ai(A + a 1/3 x 



(BIO) 



This expression is valid for x ~~ * such that a 1//3 x is not 
large, and for A/a 1 / 3 <C 1. This last condition rewrites 
L v » Ln\J Ln / L. 



(ii) The case x 



2. The limit a — ► oo : from Hermite to Airy 
function 

The first step to study this limit is to perform a rota- 
tion of +7r/4 in the complex plane of the axis of integra- 



For % — > 1 the expansion of y>(z, x) must be realized in 
the neighbourhood of — Wu/2, where the first derivative 
with respect to z vanishes (the first derivative at +^u>/2 
now diverges in the limit a — > oo). Therefore the con- 
tour of integration must be deformed in order to visit 
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the neighbourhood of z = —^Jlu/2. The new contour of 
integration is shown on the right part of figure For 
X>l/2: 



cl.r 



-i<p(x,x) — 



dz 



(Bll) 



To deal with more symmetric expressions for % > 1/2 
and x < 1/2 we remark that the contour of integration 
can also be deformed in this latter case (see the left part 
of figure[7J). For x < 1/2 : 



I 



dx 



o Vx 
Im z 



dz 

Cl+C-2 V* 



branch cut 




4 


branch cut 




— >— i 

Ci 


Rez 


-A 


C,' 














< 




1 















Im z 



(B12) 



Rez 



FIG. 7: (Color online) Left : Contour of integration in the 
complex plane for \ < 1/2 (i.e. A < 0). The vertical line 
separates the regions where Re[i(z 2 + 2Az)] is positive or neg- 
ative in the lower half of the complex plane. Right : Contour 
of integration for x > 1/2 (i.e. A > 0) 

The dominant contribution to the integral is given 
by the contribution of the segment C[. By noting that 
tp(ze~ i7r , x) = —vxu+(p{z, 1-x) for x > 1/2 and z € M + , 
we see that, for x > 1/2, 



-iv(z,x) 



-1C 



C' 



(B13) 



The solution of interest satisfies f(0) = 1 and /(l) = 0. 
As shown in the text, the harmonics of the cooperon then 
read 



C (n) (0,0) = 



1 



2V/'(0) 2 -/'(l) 2 



(C2) 



with cosh£ ff = y/jjy- 

The semiclassical approach holds for a ^> 1. In this 
case the quadratic potential acts as a high barrier in 
which the solution hardly penetrates. 

Semiclassical solution : instanton. The WKB so- 
lution of the differential equation can be expressed as 

fix) - 



Mx) 



exp± J dx'Tr(x') where the conjugate 



momentum for zero "energy" is 7r(x) = y/V(x)- Let 
us introduce : 



S(x)= [ X dx'VVte), 
Jo 



(C3) 



which is the action of the instanton penetrating inside the 
potential barrier (classical solution for imaginary time) 
for zero "energy" . We write the semiclassical solution in 
the form 



fix) 



,1/6 



R[V{ X )] 



1/4 ( 



S(x) 



+s(x) 



,(C4) 



where B sc and C sc are two coefficients to be determined 
in order to satisfy boundary conditions. 

Validity of the semiclassical approximation. The 

validity of the semiclassical approximation can be ex- 
pressed precisely by the condition Igj^Tjyl "C 1- This 
condition rewrites here as 



,1/3 



x(i-x) + 



> 1 



(C5) 



Depending on the relative magnitude of Ln and L v , there 
are two possibilities that we discuss now. We recall the 
notation A = 6/a 2 / 3 



L 2 N /L] 



Therefore, since f(x) l& dominated by C\ for x - * and 
by C[ for x - ► 1 we have 

f(x) - -ie-™f(l - x) for X -» 1 • (B14) 



APPENDIX C: SEMICLASSICAL APPROACH 

In this appendix we analyze the cooperon solution 
of (I35|l by following a semiclassical approach, valid for 
L^v <C L. We construct the solution of the equation 
[-■£2+ V(x)]f(x) = for x € [0, 1], where the potential 



V(x) = a X (l -x) + b- 



(CI) 



1. Regime L v Lm 

In this case 6/a 2 / 3 ^> 1 and the condition (|C5|) is ful- 
filled for any x G [0, 1]. We immediately find the coeffi- 
cients -B sc and C sc by imposing the boundary conditions 
and obtain : 



Therefore 



/'(I) 



f(x) = 



I njO) 8inh[5(l) - S( X )} 
tt(x) sinhS(l) 



(C6) 



/'(°) = -TT-VbcothS-^t^-Vb (C7) 



Vb 



sinh Si, 



(C8) 
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where we have introduced the notation iS; ns t = 5(1). The 
effective perimeter simply reads : 



4 ff = argcosh ~ S inst . 



2. Arbitrary L v 



(C9) 



In the case L v > Ln (i.e. b/a 2 ^ 3 < 1) the condition 
(|C5|) cannot be fulfilled near the edges of the interval. 
For x ~ the condition (|C5|) can only be satisfied for 



X > Xc whereas for % ^ 1 it is satisfied for 1 — \ > Xc- 
We have defined Xc by the two conditions a^Xc + A » 1 
and Xc "C 1 (the breakdown of the semiclassical approxi- 
mation near the edges can be simply understood by not- 
ing that, for 6 = 0, x = an d X = 1 are the turning 
points of the classical solution for imaginary time for a 
zero "energy" ) . Therefore the resolution of the differen- 
tial equation must be performed carefully near the edges. 
We separate the interval [0, 1] into three parts : 

(i) In the neighbourhood of where the potential is 
linear, the solution is a combination of two Airy func- 
tions^ : 

f(x) = B 1 Ai{aV s x + A) + Ci Bi(a x / 3 x + A) . (C10 ) 



(it) Semiclassical solution.- Sufficiently far from the 
edges of the interval, that is in the interval [\c, 1 — Xc], 
we can use the WKB solution £j4j|. 

(ra) In the neighbourhood of x — 1 the potential is al- 
most linear again and the solution reads 



f(x) = i? 2 Ai(a 1 / 3 (l-x)+A) 

+C 2 Bi(a 1 / 3 (l-x)+A). 



(Cll) 



The solution fix) is continuous and differentiable. 
Therefore the matching of the three expressions should 
now be performed in the region x ~ Xc and X ~ 1 — Xc- 
It is clear from the definition of Xc that the matching 
is realized in the Airy functions' asymptotic region. We 
obtain the following relations between the coefficients : 



(C12) 





where 



q = 2e 5,nst+ 3 A >1 



We add to these relations the conditions 

/(0) = Bi Ai(A) + C\ Bi(A) = 1 
/(l) = B 2 Ai(A) + C 2 Bi(A) = 



(C13) 



(C14) 
(C15) 



Solving these equations we find : 

q 2 Ai 



B t 



d = 



q 2 Ai 2 - Bi 2 
-Bi 



(C16) 
(C17) 



q 2 Ai 2 - Bi 2 

where Airy functions are taken at A. We eventually find : 

= /3 g 2 AiAi^-BiBi' /3 Ai^ 
W q 2 Ai 2 - Bi 2 Ai v ' 



/'(I) 



,1/3. 



qJ-K 



q 2 Ai 2 - Bi" 



,1/3. 



1 



tt q Ai' 



,(C19) 



where we have used that the Wronskian of the Airy func- 
tions is W[Ai,Bi] = AiBi' - Ai'Bi = 1/tt (see refill). 
Finally the two derivatives are 



f'iO) ~ a 
/'(I) ' 



1/3 Ai'(A) 
Ai(A) 

l/3 e -5 lnat -|A 3 / 2 



2tt Ai(A) 2 



(C20) 
(C21) 



We can check that (IC20IC21|) coincide with (|C7IC8|) in 

the limit L v <C Ijv- The effective perimeter is given by 



ff = S inst + In -47rAi(A)Ai'(A)e = 



#A 3 / 2 



(C22) 



Eq. IJC20I) corresponds to the limit l|55l) derived directly 
from the exact solution, which gives the prefactor of the 
harmonics l|59|) . 

Action of the instanton.— We analyze more into detail 
the action corresponding to the crossing of the barrier : 



5(1) = / dx^/^ 7 M 

Jo 



(C23) 



= 7TW 1 



arcsm ■ 



Vo- 



ir y/1 + a/(4b) 



(C24) 



where we recall that u> = + tA=. The action can be 

o 2 y a 

written in the form 

S inst = 5(1) = v 7 ^ 17(6/0) (C25) 
where the function rj(x) is given by : 

V(x) = I 7 + x] arctan -^L= + ^ . (C26) 
V 4 / Vix 2 

This function presents the following limiting behaviours : 
j](x) = tt Q + I ^ - I x 3/2 + 0(x 5/2 ) for x < 1 (C27) 

for x > 1 . (C28) 



^HV^+l-^ + OOr- 3 / 2 ) 
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APPENDIX D: A PERTURB ATIVE APPROACH 
TO SOLVE EQUATION (fTTTf) 

We solve equation (|37|1 with the boundary conditions 
(I41|l in the limit a<l. Let us write the solution as an 
expansion in powers of the parameter a : 



/(x) = /o(x)+/i(x)+/ 2 (x) + 



(Dl) 



where f n {x) = 0(a n )- In order to satisfy the boundary 
conditions l|41|) at any level of approximation we impose 
/ (0) = 1 and / (1) = for the order 0, and /„(0) = 
/ n (l) = for higher orders. /o(x) is solution of /q — bfo = 
0, therefore : 



/o(x) 



sinh v6(l — x) 



(D2) 



sinh yfb 

The first order term satisfies the differential equation : 

fi(x) ~ bfi(x) = a X {l ~ x)/o(x) = *(X) ■ (D3) 

The solution satisfying the appropriate boundary condi- 
tions reads : 

1 



A(x) = - w 



fo(x) I dt*(i)/ (l-t) 





+ /o(l-x) / dftf(t)/ (t) 

■'X 



(D4) 



where W 



7b 



is the Wronskian of the two solutions 



sinh \/b 

W = W[/o(x))/o(l — x)]- After some calculations the 
derivatives are found : 



f'(0) = -VbcothVb 



a 



4 sinh 2 \/b 



(D5) 



/ 1 cosh 2 \/& sinh2V6 
X ~3 + b WI^ 



+ 0{a 2 ) 



/'(I) 



-Vb 



1 



(D6) 



sinh 4 sinh 2 

1 /- coshVfo , sinh \/b \ 2 

_ coshV 5-^— +0(a 2 ) 



APPENDIX E: RELATION BETWEEN THE 
WEAK LOCALIZATION AND THE 
CONDUCTIVITY FLUCTUATIONS 

We re-examine the relation between the weak localiza- 
tion and the conductivity fluctuations studied in ref— for 
the case of the wire and used in ref—. We show that the 
relation is more general and holds for the local conductiv- 
ity a = J <J { r i *"')■ Note that it is only meaningful 
to consider a local conductivity when the distribution 
of currents is uniform (translation invariant system or a 
network with equal currents in its wires). 

The weak localization is governed only by the phase 
coherence length (Ljv and/or L v ). The study of con- 
ductivity fluctuations involves another important length 
scale : the thermal length Lt = y/D/T. Conductivity 
fluctuations are given by four contributions : the two first 
are interpreted as correlations of the diffusion constant 3 ^ 



(6a{B) 6<t(B')) (1) = 



16 /e 2 \ 2 ^D 2 r0 ° 



VoV 



3T 



dtdt' 8(t - 1') 



x drdr , r d (r,r';t)V d (r,r , ;t , y (El) 



where 5(t) = *± ( si ^ Tt ) is a function of width l/T 
and normalized to unity. The second contribution is ob- 
tained by replacing the diffuson by the cooperon V c . The 
two remaining contributions, of the form J °° dtdt' 5(t — 

t') J drdr' Re[Pd(r,r';t)Vd(r' ,r;t')], are interpreted as 
correlations of the density of states 3 ^. However, these 
two contributions are negligible^ since Lt <C Ljsr is nec- 
essary fulfilled (Lt ~ ijv corresponds to the threshold 
of strong localization). The condition Lt < i« allows 
us to replace the function S(t) by 5(t) and obtain : 



(da(B)Sa(B')) 



16 

VoP 



< j dt I drdr' 



nD 2 

V d (r,r';t) ' . (E2) 



The diffuson and cooperon are solutions of the "diffusion" 
equation 



J 



d_ 

dt 



D (V - 2ie A T f + i (Vi(r, t) - V 2 (r, t)) 



V d Jr,r';t) = 6(t)6(r-r') 



(E3) 



Azf. is the vector potential related to the magnetic field (B =p B')/2, where the sign is — for V d and + for V c . The two 
potentials V\ and V2 are the two fluctuating electric potential associated to the two conductivity bubbles. They are 
both characterized by the same fluctuations, given by the fluctuation-dissipation theorem Q17|). however Vi and V2 are 
uncorrelated since they are associated to the conductivity bubbles for two different configurations of the disorder— : 
(V,(r,t)V 3 (r',t')) v = S^TS(t - t')P d (r,r'). 

Starting from the path integral representation of the diffuson it is possible to perform the average over the fluctuating 
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potential : 



n(t)=r /-r 2 (t)=r , . , 

V Jri(0)=r' J r2 (0)=r' 



where W(r, r') was denned above by eq. (|2ip. We introduce r(r) denned for r e [0, 2t] such that r(r) = fi(r) if 
t 6 [0,t] and r(r) = r2(2t — r) if r 6 [i, 2f] . The the two path integrals can be gathered in one thanks to the 

integration over r : / dr J^ Vr x (r) f^ Q T>r 2 [T) -> /J", Q 2 * Dr(r). 



<W V d (r,r';t) 



r(2t)=r' 



Pr(r)e" 



ID Jo 



J 2t drW(r(r),r(2t-T)) 



(E5) 



y J r (0)=r' 



This leads to the following relation : 

(6a(B)5*(B'))U = -2^^(Aa (^^)> (E6) 

where the weak localization is given by I|22I23|I . Similarly 
the contribution of the cooperon gives : 

( M B) M B'))» = -2^|i(A.(^)>.(E7) 

The total correlation function (Sa(B) Sa(B')) is the sum 
of these two contributions. The relations l|E6IE7(l were 
proved by Aleiner & Blanter in ref^S by an explicit cal- 
culation of the path integral for a wire and a plane, and 
comparison to the result of AAK 4 . Here we have demon- 



strated these relations without having explicitly calcu- 
lated the path integral, which makes our proof more gen- 
eral, valid as soon as it is meaningful to consider a lo- 
cal conductivity. The important physical consequence 
of these relations is that the weak localization (i.e. the 
Altshulcr-Aronov-Spivak oscillations) and the conduc- 
tance fluctuations (i.e. the Aharonov-Bohm oscillations) 
are governed by the same length scale Ljv- For example 
we expect the amplitude of the AB oscillations to behave 
like 

/ £ \ 3/2 

6g^ B x L T ^eicp-\n\—( — ) (E8) 
for L ^> L jv ■ 
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